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Abstract. To date, over 30 multiple exoplanet systems are known, and 28% of 
stars with planets show significant evidence of a second companion. 1 briefly re- 
view these 30 systems individually, broadly grouping them into five categories: 
1) systems with 3 or more giant (Msini > 0.2Mjup) planets, 2) systems with 
, two giant planets in mean motion resonance (MMR), 3) systems with two giant 

CL(' planets not in MMR but whose dynamical evolution is affected by planet-planet 

JL ' interactions, 4) highly hierarchical systems, having two giant planets with very 

< I , large period ratios (> 30 : 1), and 5) systems of "Super-Earths", containing only 

' planets with (M sin/ < 2OM0). 

I It now appears that eccentricities are not markedly higher among planets 

in known multiple planet systems, and that planets with Msini < IMjup have 
lower eccentricities than more massive planets. The distribution of semimajor 
, axes for planets in multi-planet systems does not show the 3-day pile-up or the 

j,!^ ■ 1 AU "jump" of the apparently-single planet distribution. 

in 
a^ 

O ; 1 Introduction 
0^ 



' The fi rst multiple exoplanet system was discovered in 1992, wlien IWolszczan & Frail 



0\ '. (Il992h detected two very low-mass obj ects orbiting t he pu lsar PSR B1257H-12 using 



■ pulse timing methods. Seven years later. iButler et al. ( 1999 ) announced the first multi 



J> exoplanet system around a normal star from radial velocities, v Andromedae. 

Since then, the ever-improving precision and ever-growing temporal baseline of radial 
velocity searches has rapidly increased the number of known multiple planet systems, and 
d today at leas{3 30 such systems are knowr0. Thus, this conference comes at a special time 
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'This and the other numbers herein will be out of date at the time of this manuscript's publication - indeed, 
many were out of date by the end of the conference! At least three multiple planet systems announced shortly 
after the conference require mention, though they are not included in th e statistics here: Bou chy et al. announced 
that HD 47186 and HD 181433 harbor highly hierai'chical systems, and lMarois et al. I )2008l) imaged thi'ee young 
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in the field of multiple planet systems: there are just few enough systems now that this is 
likely the last conference where each system can still plausibly be discussed individually, 
but there are enough such that this is the first conference where we can construct some- 
thing like a statistically significant sample of planets in multiple systems to divide into 
subsamples or to compare with apparently single-planet systems. 

Multiple planet systems provide an increasingly powerful way to probe the dynamical 
origins of planets (e.g. Ipord i (i2006h ). Single-planet systems, because their orbits are 
strictly periodic, provide most of their information on the typical migration and interaction 
histories of exoplanets statistically as an ensemble. But each individual multiple planet 
system has the potential to serve as a case study of planetary system evolution. 



2 The Known Multi-Planet Systems 

I have divided the known multi-planet systems heuristically into five broad, nonexclusive 
categories: 1) systems with 3 or more giant (Msini > 0.2Mjup) planets, 2) systems with 
two giant planets in mean motion resonance, 3) systems with two giant planets not in 
MMR, but whose dynamical evolution is affected by weaker planet-planet interactions, 4) 
highly hierarchical systems, having two giant planets with very large period ratios (> 30 : 
1), and 5) systems of "Super-Earths", containing only planets with (Msini < 2OM0). The 
summaries below represent a brief and necessarily incomplete description of each system. 

2. 1 3+ Giant Planets 

Five systems are known to comprise three or more giant planets. 



2.1.1 X) Andromedae 



Butler et al. d 1999 1) announced the first multi-exoplanet system around a normal star. 



Andromedae, detected using precise radial velocity measurements from Lick Observa 
tory. The pattern of the discovery is typical of many multiplanet systems: a fter a strong, 
short-period signal from a 4.6 d, Msini = 0.68Mji,p planet was discovered dButler et al 



rt-p i 

m 



(1199m continued monitoring revealed significant structure in residuals that far exceeded 



the expected measurement uncertainties. After 2.5 y, two additional, superimposed Kep- 
lerians consistent with planets of ~ 2Mjup and 4Mjup became apparent with periods of 
- 240 d and - 1300 d. 

Interacting systems such as v Andromedae are our most powerful probes o f the dy- 
namical histories of exoplanets. To select just one example, IFord et al. I (l2005b showed 
how X) Andromedae shows good evidence that a single, strong, planet-planet scattering 
event is the origin of the modest eccentricities of the outer two planets. 



2.1.2 HD 37124 



Vogt et al. 1 (l2000h . using velocities obtained with HIRES at Keck Observatory, detected 
an apparently Jupiter-mass planet with P ~ 150 d orbiting HD 37124, a metal-poor G4 
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dwarf. Further data showe d significant deviations from the predicted velocities, and al- 
lowed iButler et al. ' ('2003') to attempt a double-planet fit including an outer planet with 
P ^6y. Gozdziewski (2003) showed that this fit was unstable, and further Keck data re- 
vealed the reason: IVogtetal.1(i2005h reported the detection of a third planet in the system, 
though with an ambiguity in its period. Today, Keck data has resolved this ambiguity. 

This system has been particularly difficult to unravel because the radial velocity am- 
plitudes of the three planets are similar and their period s are long: they have P ^ 155, 
840, and 2300 d, and Msin/ = 0.64, 0.62, and 0.68Mjup. iGozdziewski et al. I tOOdti find 
that the outermost planet's orbit is still not very well constrained, and note that in this sys- 
tem, planet-planet interactions are sufficiently strong that kinematic (sum-of-Keplerian) 
models of the radial velocities are not sufficient to describe the system. 



2.1.3 HD 74156 



Naef et al. I (120041) announced a double system with planets of Msin / — 1.86 and 6.2 Mjup 
in 52 d and 5.5 y orbits respectively around HD 74156 from data taken with ELODIE at 
Observatoire de Haute Provence. Combining published data fro m ELODIE, CORA LIE 



at La Silla Observatory, and new data acquired by HRS on HET, iBean et al. I (l2008h an 



nounced a th i rd, in termediate planet with Msin; = 0.4Mjup and a period near 1 year. 
Barnes et al. I (12008*) studied the system dynamically and found two, stable, qualitatively 
different solutions to the published RV data. They described the de tection of the third 
planet as vindication of a prediction by [Raymond & Barnes ( 2005 ) of a planet of that 
mass and orbital distance based on their "Packed Planetary System" hypothesis. 



2.1.4 ^Arae = HD 160691 



Butler et al. I (1200 Ih announced a P ^ 700 d planet with 2 Mjup orbiting /j Ara based 
on r adial velocities from UCLES on the Anglo-Australian Telescope, and soon there- 
after IJones et al. I (120021) announced that further observations revealed a linear trend in 
the resi duals, indicative of a lo ng-period, outer companion. Further AAT observations 
allowed iMcCarthv et al. I (120041) to update the fit, note that the linear trend then showed 
clear signs of curvature, and suggest a family of possible orbits including one with P ~ 
8.2 y andMsin/ — 3.1Mjiip. 



Nearly simultaneously to the latter work, ISantos et al. 1 (l2004l) . using HARPS at La 
Silla Observatory to perform high-cadence, high-precision radial velocity work, detected 
a Msin; = lO Mm compa nion in an inner 9.6 d orbit, and confirmed the outer planet. 



Finallv. iPepe et al. I ( l2007h used a combination of old and new HARPS data and the pub- 
Ushed UCLES data to detect a fourth planet with P=310 d, Msin/ = O.SMjup and deter- 
mine good orbital parameters for the outer planet for the first time in a full, dynamical, 
4-planet fit0 For the outer planet they found M sin / = 1 .8Mjup and P = 11 .5 y, and they 
revised the orbit of the b component, finding Msin / = 1 .7Mjup and P = 643 d. 



^Gozdziewski et alTI )2007l) als o announced a tentative detection of the 310-d planet from a reanalysis of 
AAT data nearly simultaneously to lPepe et al. I i20()7h . 
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Note that there is ambiguity in the literature regarding the nomenclature for these 
planets, with some authors referring to the outer planet as the e component (since it was 
the last to be characterized), and others referring to it as the c component (since it was the 
second to be detected). 



2.1.5 55 Cancri = p Cancri 



Along with Xi And b, iButler et al. I (Il997h also announced 55 Cnc h, a "51 Peg asi-type' 



planet with Msin/ = 0.84Mjup and P = 14.6 d. Further Lick data allowed Mar cy et al 
(l2002l) to announce a second planet at 5 AU with Msin/ = 4Mjup, the first extrasolar 
Jupiter analog (in terms of orbital distance). They also announced a signal from what 
appeared to be a third planet with P 45 d (roughly three times the period of the inner 
planet) and Msin / ~ 0.2Mjup, but at the time stellar rotation could not be ruled out as the 
cause. 

Using a combination of new and old Lick, EL ODIE, and HET velociti es and Hubble 



Space Telescope Fine Guidance Sensor astrometrv. lMcArthur et al. I (120041) confirmed the 
42-d planet and announced a very low amplitude 2.8-d planet with Msin/ — HM®, one 
of the first of a new class of "Hot Neptunes". Wisdom (2005), analyzing the published 
radial velocity data, challenged the reality of the 2.8 d planet, and noted a weak 260-d 
signal possibly due to a planet with Msin / ^ SOM©. 



Using a combination of new and old Lick and Keck data, iFischer et al. I (l2008h found 
a good orbital solution for all four published planets and announced the fifth, Msin/ = 
0.14Mjup, 260-d planet, making 55 Cnc the first (and, to date, only) known quintuple 
planet system. They also found that despite the near-commensurability of their orbital 
periods, the h and c components are not likely in a 3:1 mean-motion resonance because a 
dynamical integration shows that their resonant arguments do not librate. They note that, 
as with HD 37124, Keplerian models are inadequate descriptions of the existing RV data, 
and differ from the best Newtonian (dynamical) fits by > 25 m/s. 



2.2 Resonant Doubles 



These six systems contain two giant planets in or suggested to be in mean-motion reso- 
nances (MMRs). While it is difficult to understand how such planets could have formed in 
situ, differential migration could explai n how planets formed outside of resonance could 
become trapped in such an MMR (e.g. iLee and Peald (1200 II) ). The frequency and char- 
acter of such MMR systems could thus be a probe of the nature of planetary migration in 
systems with multiple giant planets. 



2.2.1 GJ 876 



Marcv et al. d 1998 1). using data from Keck Observatory, and belfosse et al.1 ( Il998l) . using 



ELODIE and CORALIE data, nearly simultaneously announced the presence of a 61-d, 
Msin / ^ 2Mjup planet orbiting GJ 876, the first know n M-dwarf planet host. After 2.5 y 
of further observations at Keck, Marcy et al. 1 (l2001al) showed that the signal was actually 
the superposition of signals from two planets in a 2: 1 mean-motion resonance, with the 
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inner planet having P= 30 d and M sin / = 0.6Mjup, making GJ 876 b and c th e first system 



clearly shown to be in a mean motion resonance. Further Keck data allowed lButler et al. 
(120041) to detect a very low-amplitude, low-mass planet in a 2.6-d orbit. 



This MMR is so strong that the orbital elements of the planets change on timescales 
shorter than the span of the extant observations of the system. The arguments of perias- 
tron of the two components precess in ^ 1 1 y, an effect clearly seen in the radial velocities 
and which complicates multi-component Keple rian fits. A dynamic al fit of the Keck ve- 
locities based on numerical integrations allowed iRivera et al. to weakly constrain 
the inclination of the system and estimate the true mass of the inner planet to be T.SM®. 



2.2.2 HD 82943 



Mavor et al. I (12004") described their discovery of the second known pair of planets in a 2: 1 
MMR orbiting HD 8294 3 (which had been publicly announced in 2000 and 2001) based 
on data from CORALIE.lLeeeLaLl (l2006l) combined the published CORALIE data with 
new Keck data to derive a dynamical solution to the system, showing that the only stable 
solutions consistent with the data are those describing a 2:1 MMR. The planets have P = 
219 and 441 d and Msin/ = 2.0 and 1.8 Mjup, respectively. 



2.2.3 HD 128311 



From an analysis of Keck data, iButler et al. I (120031) announced a P = 422 d, Msin; = 



2.2Mjup planet around the chromospherically active K dwarf HD 1 28311, for which high 
precision (< 10 m/s) can be difficult due to stellar "jitter". IVogt et al. I (l2005l) used ad- 
ditional Keck data to detect an outer companion of similar amplitude with P = 928 d 
and Msin / = 3.2Mjiip, and used dynamical simulations to show that the system is almost 
certainly locked in a 2: 1 MMR. iSandor&KlevI (l2006h suggested that since the system 
appears not to show apsidal corotation, it may owe its present state to a strong scattering 
event in its past. 



2.2.4 HD 73526 



Tinney et al. I (120031) used data from UCLES to report the detection of a P = 188 d planet 



orbiting HD 73526. Another three years of data allowed iTinney et al. I (12000) to report a 



second planet in a 378 d orbit and show t hat these planets (Msin / = 2.9 and 2.5Mjup) are 
in 2:1 resonance. ' Sandor & Kiev (2006) showed that the published solution was chaotic 
found alternative, non-chaotic (regular) orbital solutions for the system, and argued that, 
like HD 128311, the system's dynamical state showed evidence of a perturbative event 
such as a strong scattering event. 



2.2.5 HD 108874: 4:1 MMR? 



Butler et al. I (l2003h used Keck data to report a P ^ 400 d Jovian planet orbiting HD 
108874, and noted that a good fit required a linear trend be used in the model, suggesting 
a more distant companion was present in the system. By mid-2005, these residuals had 
turned over, revealing aP = 1600 d outer companion. They noted that while these orbital 
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periods are consistent with a 4: 1 MMR, that resonance is narrow and non-resonant stable 
configurations consistent with the data exist. 



2.2.6 HD 202206: 5:1 MMR? 



Udry et al. I (120021) used CORALIE data to detect a large "superplanet" or brown dwarf 
with Msin/ — 17.5Mjup orbiting HD 202206 in a 255-d orbit, on e of only a few detecti ons 



in the "brown dwarf desert". Further CORALIE data allowed ICorreia et al. I (l2005h to 
announce a second, outer companion with P = 3.8 y and Msin; = 2.4Mjup. They report 
that while the system experiences strong planet-planet interactions, stabiUty is protected 
by a 5:1 resonance. 



Th is system is especially interesting because it resembles a circumbinary planet. iGozdziewski et al. 
(1200 6') confirmed the 5:1 MMR and showed that the system is dynamically quite inter- 
esting, having different qualitative behavior depending on the inclinations of the compan- 
ions. 



2.3 Interacting Doubles 

These eight systems contain two giant planets probably not in a true mean motion reso- 
nance, but for which other planet-planet interactions can be important in the dynamical 
modeling of the system. For instance, the planets orbiting HD 12661 are in apsidal libra- 
tion. 



2.3.1 HD 12661 



Fischer et al. I (1200 Ih reported a planet orbiting the G6 dwarf HD 12661 b ased on obser- 



vatio ns at Lick and Keck Observatories. Subsequent observations allowed Fisc her et al. 



2003. ) to report that the system is actually a double. Today's best parameters ( Butler et al. 

20061)) show the b and c components having P = 263 and 1822 d, and Msin; = 2.3 and 

1.8 MjuD, respecti vely. 

fOozdziew ski & Maciejewski (2003) found the system to be near the 6:1 MMR, and 
referred to its "Janus head" of librating, anti-aligned apsidal lines. 



2.3.2 HD 155385 



Cochran et al. I (120071) used HET data to announce the lowest-metallicty planet host, the 
G subgiant HD 155385. The system is a double, with no MMR but significant planet- 
planet interactions leading to eccentricity exchange with a period of ^ 2700 y. The b and 
c components have P = 195 and 530 d, andMsin/ = 0.98 and 0.5 Mjup, respectively. 



2.3.3 HD 169830 



Naef et al. I (1200 Ih used CORALIE data to announce a P - 230 d planet with Msin ; 



2.9Mjup orbiting the F8 dwarf HD 16983C0. Further data allowed iMavor et al. I (|2004 to 



"Not to be confused with HD 69830, the triple "Super-Earth" system of j|2.5.2l 
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identif y a second, P = 5 .8 y, M sin / = 4Mjup planet in the system. iGozdziewski & Konacki 
showed that there are large exchanges of eccentricity between the components in 
coplanar configurations and that the system is stable for a large range of inclinations. 



2.3.4 HD 183263 



Marcy et al. I (l2005h used Keck data to report an Msini = 3.7Mjup plane t in a 634 d. 



eccent ric orbit around HD 183263, and note a strong residual linear trend. Wright et al. 
( 2007 ) showed that by 2007 the residuals had significant curvature, and Wright et al. 
( 120091) showed that the orbit, though still incomplete, was sufficient, in combination with 
stability analysis, to constrain the mass and period of the outer companion to P = 8.4 ± 
0.3 y and Msini = 3 — 4Mjup under the assumption that no additional companions are 
contributing to the radial velocities. 



2.3.5 47 Ursa Majoris 



The th ird star known to host exoplanets was 47 UMa, announced by iButler & Marcy 
as hosting a f = 3 y planet on a circular orbit as determined from data taken at 



Lick Observat ory. This was t he firs t planet strongly reminiscent of the gas giants in our 



Solar System. iFischer et al. I (l2002h studied another 5 years of Lick data and reported a 

sec ond planet in the sy s tem, h aving P 1 100 d and Msin/ = 2.5Mjup. Naef et al. 

(200i 



and Iwittenmyer et al. I (2007), using ELODIE and HET data, respectively, have noted 
that a second planet with the reported parameters is not apparent in their data. Future 
observations should clarify the situation. 



2.3.6 HIP 14810 



Butler et al. I (12006b included HIP 14810 b in the Catalog of Nearby Exopl anets based on 
prelim inary data from Keck Observatory taken as part of the N2K survey d pischer et al. 



(120051) 1 and further observations allowed Wright et al. (2007) to provide an orbital solu- 
tion for two planets, having P = 6.7 and 95 d and Msini = 3.9 and 0.76 Mjup. Because 
the outer planet's orbit is currently poorly sampled, further observations will help clarify 
the nature of this system and allow for it to be better studied dynamically. 



2.3.7 OGLE-2006-BLG-109L 



Gaudi et al. I (120081) announced the remarkable detection of a double-planet system around 
the distant (d = 1.5kpc) star OGLE-2006-BLG-109L during a microlensing event. Due 
to the great distance to this system and the non-repeating nature of microlensing de- 
tections, the study of the d ynamics of this system with specificity is difficult (but see 



Malhotra & Minton I ( 120081) ). This detection demonstrates the promise of microlensing 
as a method to build up statistics of multi-planet systems, including those composed of 
rocky planets at a few AU. The most likely masses and orbital distances of these planets 
are m ^ 0.71 and 0.27Mjup and a ^ 2.3 and 4.6 AU, making this system around a K star 
a "scaled-down Juipter-Saturn analog" 
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2.3.8 HD 102272 



Niedzielski et al. I (120081) used HET to announce two companions orbiting the giant star 
HD 102272. The inner, Msin/ — 5.9Mjup planet orbiting at 0.6 AU is the closest known 
companion to a star with M > 1.5Mq. Correlated residuals to a one-planet fit indicate a 
second planet of uncertain orbital period and mass, but the velocities are consistent with 
a period of P ^ 520d. 



2.4 Highly Hierarchical Doubles 

These seven systems have orbital period ratios greater than 30; 1, and so have little interac- 
tion between their components and can usually be well-modeled without resort to A^-body 
simulations. 



2.4.1 HD 168443 



From Keck data, iMarcy et al. 1 (1200 ih discovered a pair of massive, highly hierarchical 
companions orbiting HD 168443: a close-in P = 58 d innerplanet with Msin/ = 8. 16Mjup 
and an outer object with P = 4.8 y and Msin; = 18.4Mjup. Like HD 202206, this sys- 
tem has wide separation between its components and contains a "super-planet", but here 
the orbit of the lighter obje ct is of the"S-type" (it orbits one, not both of the massive 
companions. iDvorak I (Il983h '). 



2.4.2 HD 187123 



Using Keck data. lButler et al. I (Il998h announced an M sin ; = O.SMju n planet in a 3-dav or- 
bit around HD 187123, a solar "twin". Further observations allowed lWright et al. I ( 120071) 
to report the existence of an outer companion with orbital period > 10 y, and to constrain 
its minimum m ass to be plane t ary. U sing subsequent observations of the apparent clos- 
ing of the orbit, Wright et al. ( 20091) constrained its orbit to have P = 10.5±0.5 y and 
Msin/ ~ 2.0±0.1Mjup, under the assumption that no other planets are influencing the 
observations. 



2.4.3 HD 68988 



[Voit et alj ( 120021) used Keck data to report a 6. 3-d, Msin / ~ 2Mjup planet orbiting HD 



68988. Subsequent observ ations have revealed a long-period outer companion of uncer- 
tain mass and period, with IWright et al. I (l2007h constrained to be 1 ly < P < 60 y and 
6Mjup < Msin / < 20Mjup. 



2.4.4 HD 38529 



Fischer et al. I ( 1200 Ih used Lick and Keck data to report a 14. 3-d, Msin/ = 0.8Mjup com- 



panion to the old G subgiant HD 38529, and not ed that the residuals t o the fit suggested an 



outer companion. Further observations allowed iFischer et al. 
P = 5.9 y, Msin/ 



(l2003h to confirm an outer. 



13.2Mjup companion. Spitzer observations by iMoro-Martm et 
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(l2007b reveal that HD 38529 has infrared excess consistent with dust-producing planetes- 
simals at around 5 AU (exterior to both planets). 



2.4.5 HD 217107 

Fischer et al. 1 (ll999h used Lick data to discover the first, inner P 



planet around the G7 dwaif HD 217107. IVogt et al 



7.1 d,Msin/= 1.3Mjup 
(l2005h used new and old Lick and 



Keck data to identify a long per iod outer companion of uncertain mass and period. Further 
observations at Keck allowed IWright et al. [ (120091) to constrain the outer companion's 
mass and period to be f ^ 11.7 y and Msin/ ^ 2.6Mjup under the assumption that no 
other planets are contributing to the radial velocities. 



2.4.6 HD 11964 



Butler et al. I (l2006l) announced a M sin ; — 0.6Mjup planet orbiting the slightly evolved star 



HD 1 1964 in a 5.3 y orbit. Subsequent observations suggested a l ow mass inner plane t 



with P = 38 d (IWright et al. I ( 2007) ). and further monitoring allowed lWright et al. 
to confirm this signal as being to due a Msin / = 23M0 companion. 



(2009) 



2.4.7 GJ 777 A = HD 190360 

Naef et al. I (l2003b used data from ELODIE and the AFOE spectrograph at Whipple Ob- 
servatory to detecil a long-period, Jovian planet around GJ 777 A. .Vogt et al. ( 2005) used 
Keck data to confirm the outer planet and revise its orbital parameters, and to announce a 
second, close-in, low-mass companion. The b and c components have P = 8.0 y and 17.1 
d, and Msin; = 1.6Mjup and 19M©. 



2.5 Systems of "Rocky" Planets and "Super-Earths" 

These four systems compose the present-day bookends of multiple-planet systems. The 
pulsar triple-planet system is the first exoplanetary system known and remains a fascinat- 
ing example of the limits of multiple-planet detection. The latest frontier of exoplanet 
research with radial velocities is the hunt for rocky planets, and these latest detections 
of multiple "Super-Earths", all from the HARPS spectrograph, represent the penultimate 
step toward the definitive detection of rocky planets. 

Because they do not transit, the actual compositions and masses of these planets is 
unknown, and so the monicker "rocky" is probably only truly appropriate for the pulsar 
planets. 



2.5.1 PSRB1257H-12 



Wolszczan & Frail I (119921) detected two very low-mass (Msin; — 3.9 and 4.3 M^) ob- 
jects orbiting the pulsar PSR B1257H-12 using pulse timing methods. Two years later. 



'a preliminary orbit for this system also appears in lUdrv et al. I J2003h . 
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Wolszczan I (1 1 994l) detected the planet-planet interactions from the planets' 3:2 MMR, as 
well as a smaller, 0.02 object in the system. 



2.5.2 HD 69830 



Based on observations with HARPS, iLovis et al. I (120061) announced a "triple-Neptune" 
system orbiting the K dwarf HD 69830, having periods o f ^ 9, 32, and 200 days. This 



system is especially interesting because iBeichman et al. I (l2005h had used Spitzer pho- 
tometry and spectroscopy to reveal an infrared excess characteristic of a large cloud of 
fine silicate dust within a few AU of the star, suggestive of a large asteroid belt or "super- 
comet". 



2.5.3 GJ581 



Also using HARPS, iBonfils et al. I (l2005h announced a Msini = 0.052Mjup planet in a 



5.366 d orbit around GJ 581, making that star only t he third M dwarf k nown to host a 



planet. After collecting additional data, two years later lUdrv et al. I (l2007h announced the 



discovery of two additional planets in the system with minimum masses of 5 and 7.7 M®, 
with pe riods of 12.9 and 8 3.6 d. The outermost planet in the system sits at 0.25 AU, a 
location lUdrv et al. I (120071) a ndlSe lsis et al. (2007) identify as being near the "cold edge" 

dKasti 



of the star's Habitable Zone dKasting et al. I (tl993|)) 



2.5.4 HD 40307 



Mavor et al. I (l2009h announced a third triple "Super-Earth" system from HARPS around 
the K2, dwarf HD 40307, and noted a linear trend in the radial velocities suggestive of 
a fourth, outer companion, as well. The low metallicity of this star ([Fe/H] = —0.31) 
and others has led to the suggestion that the well known metallicity dependence of planet 
occurrence breaks down for low-mass stars. 



2.6 Future Detections and the Multiplicty Rate 

Many systems show strong evidence of second planets due to long-period companions 
whose orbits are too incomplete for strong constraints to be put on their masses. For 
instance, 14 Her c l early has a long period companion of some so rt, p robably planetary 
(iNaef et al. I(l2004 . lwitt enmver et al.1(l2007l).IWright et al. I(l2007b . and lOozdziewski et al. 
(boos!)), as may GJ 317 (Johnso n et al. I (120071) ). 



Of the 200 planet-bearing stars within 200 pc, the 28 above (not counting PSR B 1 257h- 1 2 
and OGLE-2006-BLG-109L) constitute 14% of the total. An additional 27 (including 14 
Her and GJ 3 17) show clear evidence of trends in their residuals and no evidence of stellar 
duplicity, meaning that the true multiplicity rate is at least 28%. 



3 Statistical Properties of tlie l\/lultiple Planet Sample 



Wright et al. (12009 ^ contains a catalog of the above multi-planet systems including up- 



dated latest orbital parameters for 10 of the systems. From this catalog, several intriguing 
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patterns emerge; these features represent an opportunity to test models of planet forma- 
tion, migration, and the origin of eccentricities: 

• Planets in multiple-planet systems have eccentricities no higher than single planets 
(see Figure[T]). 

• The distribution of orbital distances of planets in multi-planet systems and single 
planets are inconsistent: single-planet systems show a pile-up at f ^ 3 days and a 
jump near 1 AU, while multi-planet systems show a more uniform distribution in 
log-period (see Figure |2]i. 



30 

25 

20 

15 

10 

5 




Single-planets 
Multi-planets 



0.0 0.2 0.4 0.6 

e 



0.8 



1.0 



Fig. 1. Distribution of eccentricities of exoplanets for known multiple planet systems (solid) and 
apparently single planet systems (dashed). Note the high eccentricity orbits, e > 0.6 occur predom- 
inantly in single planets. 



In addition, among all planetary systems: 

• There may be an emerging, pos itive correlation betw een stellar mass and giant- 
planet semi-major axis (see, e.g. Jjohnson et al. I OOOTD V 



• Exoplanets more massive than Jupiter have eccentricities broadly distributed across 
< e < 0.5, while lower-mass exoplanets exhibit a distribution peaked near e = 
(see Figure Is). 
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Fig. 2. Distribution of semimajor axes of exoplanets for multiple-planet systems (solid) and ap- 
parently single systems (dashed). Note the enhanced frequency of hot jupiters and the jump in 
abundance beyond 1 AU in the single-planet systems. 



From Figure |2] it is clear that the orbital distances of planets in systems currently 
known to be multiple are not drawn from the apparently single-planet distribution, indi- 
cating that the migration mechanisms operate differently in these populations. Yet the 
eccentricity distributions of these two populations are nearly identical (Figure [U, sug- 
gesting that the mechanisms of eccentricity excitation are similar. Figure |3] shows that 
whatever those mechanisms are, they are ultimately more effective in pumping the eccen- 
tricities of planets with M sin/ > IMjup. Reproducing these distributions will be a good 
test for future models and simulations of planet formation and dynamical evolution. 



References 

Barnes, R., Gozdziewski, K., & Raymond, S. N. 2008, Astrophysical Journal Letters, 
680, L57 

Bean, J. L., McArthur, B. E., Benedict, G. F, & Armstrong, A. 2008, Asti-ophysical 
Journal, 672, 1202 



A Survey of Multiple Planet Systems 



13 



20 



15 



10 - 



M > 1.0 M 
M < 1.0 M 



Jup 
Jup 



5- 







J 



I 



1 



0.0 0.2 0.4 0.6 

Eccentricity 



0.8 



1.0 



Fig. 3. Distribution of eccentricities of exoplanets with M sin i < l.OMjup (solid) Msini > l.OMjyp 
(dashied). The tidally circularized hot jupiters have been removed. Note that the eccentricity of 
planets of minimum mass < 1 .OMjup peaks at eccentricity < 0.2, while the eccentricities e of planets 
of minimum mass > l.OMjup are distributed broadly from 0.0 < e < 0.5. 



Beichman, C. A., Bryden, G., Gautier, T. N., Stapelfeldt, K. R., Werner, M. W., Misselt, 
K., Rieke, G., Stansberry, J., Trilling, D. 2005, Astrophysical Journal 626, 1061-1069. 

Bonfils, X., Forveille, T., Delfosse, X., Udry, S., Mayor, M., Perrier, C, Bouchy, R, Pepe, 
F., Queloz, D., Bertaux, J.-L. 2005, Astronomy and Astrophysics 443, L15-L18. 

Butler, R. P., & Mai'cy, G. W. 1996, Astrophysical Journal Letters, 464, L153 

Butler, R. P, Mai'cy, G. W., Williams, E., Hauser, H., Shirts, P 1997. Three New "51 
Pegasi-Type" Planets. Astrophysical Journal 474, LI 15. 

Butler, R. P., Marcy, G. W., Vogt, S. S., & Apps, K. 1998, Publications of the Astronomi- 
cal Society of the Pacific, 110, 1389 

Butler, R. P, Marcy, G. W., Fischer, D. A., Brown, T. M., Contos, A. R., Korzennik, S. G., 
Nisenson, P, Noyes, R. W. 1999, Astrophysical Journal 526, 916-927. 



14 



Title : will be set by the publisher 



Butler, R. P., Tinney, C. G., Marcy, G. W., Jones, H. R. A., Penny, A. J., & Apps, K. 2001, 
Astrophysical Journal, 555, 410 

Butler, R. P, Marcy, G. W., Vogt, S. S., Fischer, D. A., Henry, G. W., LaughUn, G., & 
Wright, J. T. 2003, Astrophysical Journal, 582, 455 

Butler, R. P, Vogt, S. S., Marcy, G. W., Fischer, D. A., Wright, J. T., Henry, G. W., 
Laughlin, G., & Lissauer, J. J. 2004, Astrophysical Journal, 617, 580 

Butler, R. P, et al. 2006, Astrophysical Journal, 646, 505 

Cochran, W. D., Endl, M., Wittenmyer, R. A., & Bean, J. L. 2007, Astrophysical Journal, 
665, 1407 

Correia, A. C. M., Udry, S., Mayor, M., Laskar, J., Naef, D., Pepe, F., Queloz, D., & 
Santos, N. C. 2005, Astronomy & Astrophysics, 440, 751 

Delfosse, X., Forveille, T., Mayor, M., Perrier, C, Naef, D., & Queloz, D. 1998, Asti'on- 
omy & Astrophysics, 338, L67 

Dvorak, R., 1983, Osterreichische Akademie der Wissenschaften, Mathematisch- 
Naturwissenschaftliche Klasse, Bd. 191, p. 423 

Fischer, D. A., Marcy, G. W., Butler, R. P, Vogt, S. S., & Apps, K. 1999, Publications of 
the Astronomical Society of the Pacific, 1 1 1, 50 

Fischer, D. A., Marcy, G. W., Butler, R. P, Vogt, S. S., Frink, S., & Apps, K. 2001, 
Astrophysical Journal, 551, 1107 

Fischer, D. A., Marcy, G. W., Butler, R. P, Laughlin, G., & Vogt, S. S. 2002, Astrophysi- 
cal Journal, 564, 1028 

Fischer, D. A., Mai'cy, G. W., Butler, R. P, Vogt, S. S., Henry, G. W., Pourbaix, D., Walp, 
B., Misch, A. A., Wright, J. T. 2003, Astrophysical Journal 586, 1394-1408. 

Fischer, D. A., et al. 2005, Astrophysical Journal, 620, 481 

Fischer, D. A., et al. 2008, Astrophysical Journal, 675, 790 

Ford, E. B., Lystad, V, Rasio, F A. 2005, Nature 434, 873-876. 

Ford, E. B. 2006, New Horizons in Astronomy: Frank N. Bash Symposium 352, 15. 

Gaudi, B. S., et al. 2008, Science, 319, 927 

Gozdziewski, K., & Maciejewski, A. J. 2003, Astrophysical Journal Letters, 586, L153 
Gozdziewski, K. 2003, Astronomy & Astrophysics, 398, 315 
Gozdziewski, K., & Konacki, M. 2004, Astrophysical Journal, 610, 1093 



A Survey of Multiple Planet Systems 



15 



Gozdziewski, K., Konacki, M., & Maciejewski, A. J. 2006, Astrophysical Journal, 645, 
688 

Gozdziewski, K., Maciejewski, A. J., Migaszewski, C. 2007. On the Extrasolar Multi- 
planet System around HD 160691. Astrophysical Journal 657, 546-558. 

Gozdziewski, K., Migaszewski, C, & Konacki, M. 2008, Monthly Notices of the Royal 
Astronomical Society, 385, 957 

Johnson, J. A., Butler, R. R, Mai'cy, G. W., Fischer, D. A., Vogt, S. S., Wright, J. T., & 
Peek, K. M. G. 2007, Astrophysical Journal, 670, 833 

Jones, H. R. A., Paul Butler, R., Marcy, G. W., Tinney, C. G., Penny, A. J., McCarthy, C, 
& Carter, B. D. 2002, Monthly Notices of the Royal Astronomical Society, 337, 1 170 

Kasting, J. R, Whitmire, D. P, & Reynolds, R. T. 1993, Icarus, 101, 108 

Lee, M. H., Peale, S. J. 2001. Evolution of the GJ876 Planets into the 2:1 Orbital Reso- 
nance. Bulletin of the American Astronomical Society 33, 1 198. 

Lee, M. H., Buder, R. P, Fischer, D. A., Marcy, G. W., & Vogt, S. S. 2006, Astrophysical 
Journal, 641, 1178 

Lovis, C, et al. 2006, Nature, 441, 305 

Malhotra, R., & Minton, D. A. 2008, Astrophysical Journal Letters, 683, L67 

Marcy, G. W., Butler, R. P, Vogt, S. S., Fischer, D., & Lissauer, J. J. 1998, Astrophysical 
Journal Letters, 505, L147 

Marcy, G. W., Buder, R. P, Fischer, D., Vogt, S. S., Lissauer, J. J., & Rivera, E. J. 2001, 
Astrophysical Journal, 556, 296 

Marcy, G. W., Buder, R. P, Vogt, S. S., Liu, M. C, LaughHn, G., Apps, K., Graham, J. R., 
Lloyd, J., Luhman, K. L., Jayawardhana, R. 2001, Astrophysical Journal 555, 418-425. 

Marcy, G. W., Buder, R. P, Fischer, D. A., Laughlin, G., Vogt, S. S., Henry, G. W., & 
Pourbaix, D. 2002, Astrophysical Journal, 581, 1375 

Marcy, G. W., Butler, R. P, Vogt, S. S., Fischer, D. A., Henry, G. W., LaughUn, G., 
Wright, J. T., & Johnson, J. A. 2005, Astrophysical Journal, 619, 570 

Mayor, M., Udry, S., Naef, D., Pepe, F, Queloz, D., Santos, N. C, & Burnet, M. 2004, 
Astronomy & Astrophysics, 415, 391 

Marois, C, Macintosh, B., Barman, T., Zuckerman, B., Song, I., Patience, J., Lafreniere, 
D., Doyon, R. 2008, Science 322, 1348. 

Mayor, M., Udry, S., Lovis, C, Pepe, F, Queloz, D., Benz, W., Bertaux, J.-L., Bouchy, 
F., Mordasini, C, Segransan, D. 2009, Astronomy and Astrophysics 493, 639-644. 



16 



Title : will be set by the publisher 



McArthur, B. E., et al. 2004, AstrophysicalJournal Letters, 614, L81 

McCarthy, C, Butler, R. P., Tinney, C. G., Jones, H. R. A., Marcy, G. W., Carter, B., 
Penny, A. J., & Fischer, D. A. 2004, Astrophysical Journal, 617, 575 

Moro-Martm, A., Malhotra, R., Carpenter, J. M., Hillenbrand, L. A., Wolf, S., Meyer, 
M. R., Hollenbach, D., Najita, J., Henning, T. 2007, Astrophysical Journal 668, 1165- 
1173. 

Naef, D., Mayor, M., Pepe, R, Queloz, D., Santos, N. C, Udry, S., & Burnet, M. 2001, 
Astronomy & Astrophysics, 375, 205 

Naef, D., et al. 2003, Astronomy & Astrophysics, 410, 1051 

Naef, D., Mayor, M., Beuzit, J. L., Perrier, C, Queloz, D., Sivan, J. P, & Udry, S. 2004, 
Astronomy & Astrophysics, 414, 351 

Niedzielski, A., Gozdziewski, K., Wolszczan, A., Konacki, M., Nowak, G., & Zielinski, 
P 2008, AstrophysicalJournal, 693, 276 

Pepe, E, et al. 2007, Astronomy & Astrophysics, 462, 769 

Raymond, S. N., & Barnes, R. 2005, Astrophysical Journal, 619, 549 

Rivera, E. J., et al. 2005, Astrophysical Journal, 634, 625 

Sandor, Z., & Kley, W. 2006, Astronomy & Asti'ophysics, 451, L31 

Sandor, Z., Kley, W., & Klagyivik, P 2007, Astronomy & Astrophysics, 472, 981 

Santos, N. C, et al. 2004, Astronomy & Astrophysics, 426, L19 

Selsis, P., Kasting, J. P., Levrard, B., Paillet, J., Ribas, 1., & Delfosse, X. 2007, Astronomy 
& Astrophysics, 476, 1373 

Tinney, C. G., Butler, R. P, Marcy, G. W., Jones, H. R. A., Penny, A. J., McCarthy, C, 
Carter, B. D., & Bond, J. 2003, Astrophysical Journal, 587, 423 

Tinney, C. G., Butler, R. P, Marcy, G. W., Jones, H. R. A., LaughUn, G., Carter, B. D., 
Bailey, J. A., & O'Toole, S. 2006, AstrophysicalJournal, 647, 594 

Udry, S., Mayor, M., Naef, D., Pepe, E, Queloz, D., Santos, N. C, & Burnet, M. 2002, 
Astronomy & Astrophysics, 390, 267 

Udry, S., Mayor, M., & Queloz, D. 2003, Scientific Frontiers in Research on Extrasolar 
Planets, 294, 17 

Udry, S., and 10 colleagues 2007, 469, L43-L47. 

Vogt, S. S., Marcy, G. W., Butler, R. R, & Apps, K. 2000, AstrophysicalJournal, 536, 902 



A Survey of Multiple Planet Systems 



17 



Vogt, S. S., Butler, R. P., Marcy, G. W., Fischer, D. A., Pourbaix, D., Apps, K., & Laugh- 
lin, G. 2002, Astrophysical Journal, 568, 352 

Vogt, S. S., Butler, R. R, Marcy, G. W., Fischer, D. A., Henry, G. W., LaughHn, G., 
Wright, J. T., & Johnson, J. A. 2005, Astrophysical Journal, 632, 638 

Wisdom, J. 2005, Bulletin of the American Astronomical Society, 37, 525 

Wittenmyer, R. A., Endl, M., & Cochran, W. D. 2007, Astrophysical Journal, 654, 625 

Wolszczan, A., Frail, D. A. 1992, Nature 355, 145-147 

Wolszczan, A. 1994, Science 264, 538 

Wright, J. T., and 10 colleagues 2007, Astrophysical Journal 657, 533-545. 

Wright, J. T., Upadhyay, S., Marcy, G. W., Fischer, D. A., Ford, E. B., Johnson, J. A. 
2009, Astrophysical Journal 693, 1084-1099. 



